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We propose a method for analyzing two-dimensional symmetry protected topological (SPT) wave- 
functions using a correspondence with conformal field theories (CFTs) and integrable lattice models. 
This method generalizes the CFT approach for the fractional quantum Hall effect wherein the wave- 
function amplitude is written as a many-operator correlator in the CFT. Adopting a bottom-up 
approach, we start from various known microscopic wavefunctions of SPTs with discrete symme¬ 
tries and show how the CFT description emerges at large scale, thereby revealing a deep connection 
between group cocyles and critical, sometimes integrable, models. We show that the CFT describing 
the bulk wavefunction is often also the one describing the entanglement spectrum, but not always. 
Using a plasma analogy, we also prove the existence of hidden quasi-long-range order for a large 
class of SPTs. Finally, we show how response to symmetry fluxes is easily described in terms of the 
CFT. 


I. INTRODUCTION 

In the past three decades, topological phases of mat¬ 
ter have attracted a large amount of interest due to their 
tendency to exhibit highly robust quantum phenomena 
which have various applications in quantum engineering 
and metrology. One of the current frontiers in the field 
aims at understanding the variety of novel topological 
phases which arise when some extra symmetries are not 
allowed to be broken. For example, it was shown m 
that, given an Ising {Z 2 ) symmetry, there are two topo¬ 
logically distinct Ising paramagnets in two dimensions. 
This can be thought of as the spin analog of topological 
insulators [4] and accordingly the topological Ising para- 
magnet must have gapless magnetic excitations on its 
boundary. While numerous types of topological insula¬ 
tors have been realized experimentally HU, such “topo¬ 
logical paramagnets” have thus far only been realized in 
the ID context [7], although several ideas have been put 
forward concerning 2D H [5] and 3D m- 

Since this discovery, a variety of topological param¬ 
agnets (more commonly known as short-range entan¬ 
gled bosonic SPTs) with different symmetries and dif¬ 
ferent dimensions have been explored using various ad¬ 
vanced tools HEHIS]. Long-range entangled versions 
of these phases (coined symmetry enriched topological 
(SET) phases) have been studied as well [T6l[^. Still, 
open questions remain concerning the scope of these ap¬ 
proaches their extensions to fermionic systems |18j . 
and their relevance to experimentally feasible models. 
Also, microscopic lattice models realizing fractional topo¬ 
logical paramagnets [71] (for instance gapped spin sys¬ 
tems with an Ising {Z 2 ) symmetry supporting bulk exci¬ 
tations with a Z 4 symmetry) are very scarce [75] . 

A powerful and conceptually simple tool in analyzing 
the fractional quantum Hall effect (FQHE), a prominent 
2D topological phase, is the FQHE-CFT correspondence 
[I9ll2^ . While CFTs are usually used to describe critical 
states, here the ground state wavefunction of a 2D topo¬ 
logical phase is expressed as a many-operator correlation 


function in the holomorphic (chiral) part of a CFT: 

Ii/'fqhe) = X ( (1) 

where Zi are the particle positions given on the complex 
plane and the expectation value is computed in a given 
CFT in which is a given operator. Choosing differ¬ 
ent CFTs and different operators O yields a variety of 
Abelian, non-Abelian, bosonic, and fermionic fractional 
quantum Hall (FQH) phases. The quasiparticle statis¬ 
tics, ground state degeneracies, and edge spectrum are 
all readily deduced from known properties of CFTs [22] . 
Furthermore microscopic Hamiltonians which stabilize 
these phases can be written down. 

There are many recent works showing similarities be¬ 
tween SPTs and the FQHE [Tn[23H25] . Most relevant to 
our discussion is Ref. m where SPT phases are analyzed 
through the prism of bi-layer or multi-layer quantum Hall 
effects. Using the K-matrix formulation, several classes 
(though not all [75]) of SPTs were shown to correspond to 
compact boson CFTs with an unconventional (chiral) ac¬ 
tion of the physical symmetry on the left and right goers 
(see also Ref. [M])- While this suggests that some ground 
state wavefunctions of SPT phases can be presented as 
correlators in these compact boson CFTs, this continuum 
observation appears physically relevant only for bi-layer 
quantum Hall setups and not for spins on a lattice. The 
main microscopic approach to studying SPTs on a lat¬ 
tice is the group cohomology approach]!] [5] and, at least 
for bosonic SPTs, it is also more comprehensive in its 
scope. Physically, however, it is somewhat opaque and 
it is in particular unclear what CFTs one can associate 
with these wavefunctions and whether this is beneficial 
in some way. 

Our work therefore provides a distinct and unify¬ 
ing approach to studying SPT phases which general¬ 
izes the above CFT-FQHE correspondence and at the 
same time interpolates between the continuum field the¬ 
ory approaches of Refs. m [l5] and the microscopic 
group cohomology approach. It is shown that micro- 
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scopic 2D bosonic SPT wavefunctions, when written in 
the symmetry-charge basis (i.e. the basis on which the 
symmetry acts diagonally), appear as a many-operator 
correlation function in a non-chiral CFT (i.e. where 
both holomorphic and anti-holomorphic parts are taken). 
Based on loop models and notions of discrete flux at¬ 
tachment |25j . a microscopic toolbox is developed which 
allows us to identify CFTs associated with a large va¬ 
riety of SPTs including SPTs for which the continuum 
approaches do not apply in any straightforward way. Fur¬ 
thermore, the expected symmetry-flux responses [521 HI] 
are obtained. For a Zjv >2 symmetry, this CFT ap¬ 
proach is used to establish the presence of hidden or¬ 
der |25l|27l|28], a property which to the best of our knowl¬ 
edge has not been derived before. Somewhat surprisingly, 
unlike in the CFT-FQHE context where it is believed to 
be impossible to obtain an exact tensor product state 
(TPS) description, here there appears to be no tension 
between having a TPS and writing the wave function as 
correlators in a CFT. Lastly, we show that for all group 
cohomology wave functions, the entanglement spectrum 
is given by the spectrum of the CFT. 


II. MICROSCOPIC SPT-CFT 
CORRESPONDENCE 


We consider the ground state \'ip) of a 2D bosonic SPT 
based on a discrete symmetry group G. The Hilbert 
space is given by degrees of freedom G G lying on the 
sites of a triangular lattice at position r. For simplicity, 
let us focus on G = Z^, where we label group elements 
as the numbers {0,..., iV— 1} ~ G. In the basis Ipr), the 
action of g G G is given by g\(j)r) = \g + (pr)- Since it ro¬ 
tates under the action of G, we refer to this basis as the 
symmetry-phase basis. The SPT wavefunction is then 
given by l'^) = |{<?^r}) where all the 

Af configurations {4>r} Etre summed over. For all known 
wavefunctions, is given by a product of local phase 

factors [21 [55]. 

To establish similarities with FQH wavefunctions, 
which are written as a function of the charge positions 
(the charge being electrons in that case), it is natural to 
use the symmetry-charge basis, G {0,..., iV — 1}, on 
which g acts diagonally: g\ar) = [5S]. For 

non-Abelian G, this would be the representation basis. 
It will be advantageous to enumerate the {ur} basis by 
a set of unordered positions {±i, r^} of “charge ±1 parti¬ 
cles” such that ar = mod N. Redundancy in 

this description is removed by always taking the smallest 


number of particles. The SPT wavefunction reads 

{At} i 

MAt} 

{At} 

The amplitudes of \'ip) in the symmetry-charge basis ap¬ 
pear as correlators in a statistical mechanics model whose 
configuration space is {4>r}^ whose “Boltzmann” weights 
are and where an operator 0 ±{r) = gG'^'^iAT/N 

inserted at every charge ±1 particle position. 

A major motivation for our approach is the observation 
that two-point correlation functions within this partition 
function, also known as strange correlators between the 
SPT state and an ideal paramagnet, appear to be criti¬ 
cal [30] . In the following, we will show explicitly for many 
examples that this Z theory is indeed critical, and that 
it also gives the entanglement spectrum theory. Even so, 
for Eq. [^ to be a useful generalization of Eq. [^ three 
main issues must be addressed: (a) the identification of 
the CFT governing this critical behavior, (b) the associa¬ 
tion between 0± and CFT operators, and (c) the strong 
ultraviolet corrections to CFT predictions, coming from 
the fact that the 0 ± operator insertions are dense on the 
lattice scale. 

To address the last issue, it will be advantageous to 
consider the dilute symmetry charge regime where the 
number of particles per lattice site p is made small by 
adding a particle fugacity factor to the wavefunction. 
In this regime, the correlation function involves 

mainly long-range features which can therefore be stud¬ 
ied using CFT. Provided that the wavefunction near 
p = 0 is adiabatically connected to that with p « 1, the 
qualitative features thus obtained should persist down 
to the lattice scale (p « 1). To establish this adiabatic 
connection, in Appendix [C] we use the fact that \il}) can 
be written as a tensor product state (TPS) to show that 
a continuous family of local fugacity-dependent Hamil¬ 
tonians exists for which these fugacity-dependent wave- 
functions are the exact ground states. Moreover, these 
ground states are found to be strictly short-range corre¬ 
lated, strongly implying that this family of Hamiltonians 
is gapped as required. Note that this adiabatic connec¬ 
tion is akin to that between FQH phases in the continuum 
(p —>■ 0) and their lattice counterparts, coined Fractional 
Chern Insulators, for which p is of order one. Here again, 
there is ample numerical [211 - 124] and analytical [25U29] 
evidence. 


III. ENTANGLEMENT SPECTRUM 

We now explicitly show the relation between the the¬ 
ory describing the SPT bulk wavefunction (Z) and the 
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FIG. 1: Bipartition of the infinite cylinder (infinite in the x 
direction, periodic boundary conditions along y with width 
L). <E>± is an index running over all possible conhgurations of 
<(>( 0 =^, 2 /). 


theory describing the entanglement spectrum. While it 
is difficult to extract this spectrum in general, a signif¬ 
icant simplification arises from an important feature of 
all the “fixed-point wave functions” [53] we consider in 
this work: ^{ 0 ,,} always appears as a product over tri¬ 
angles of local phase factors. The partition function for 
a triangular lattice on a torus of length N and width L 
(see figure can therefore be written as 

Z = Tt[{TiT2 )^/% (3) 

where Ti and T 2 are \G\^ x \G\^ transfer matrices for 
which the matrix elements are given by the prod¬ 

uct of phase factors for a given column of triangles and 
where the group elements on the left and right side of 
this column are given by <i> and respectively. If Z de¬ 
scribes a conformal field theory, the low lying spectrum 
of log(TiT 2 ) is, up to a simple normalization, equal to 
the spectrum of scaling dimensions in the CFT. 

Let us consider a system defined on an infinite cylinder 
and consider cutting the cylinder in half. We define r~ S 
{{x,y)\x < 0"}, r+ e {lx,y)\x > 0+}, (j)± G 
G}.We compute the reduced density matrix of the left 
half of the cylinder. It is given by 

^ Ip Ip* [(/)'_, (1)+] (4) 

<P + 


leading to 

<s>+ 

= A_[^_]A*_[^'_]%,_^^,_ ( 6 ) 

= V'A(«!>-)V’A(</>-) 

A 

where T = TT\ 7$,$' u^, = X and ?(’>((/)_) = 
A-[4>J\u^_. Using the fact that |A_[(/>_]p = 1, it is 
easy to see that the ipx form an orthonormal set. We can 
therefore conclude that the entanglement spectrum A is 
given by the spectrum of T. 

We now have two theories defined by two different 
transfer matrices: the theory for the bulk wavefunction 
is given by T 1 T 2 and the theory for the entanglement 
spectrum is given by TiT\. Now, within the group co¬ 
homology construction of SPT wavefunctions, the local 
cocycle factor is taken with a complex conjugate for, say, 
all left-pointing triangles. It is easy to see that this com¬ 
plex conjugation ensures that T 2 =t\, and that the two 
theories are therefore the same. 

While the bulk and edge theory are the same for group 
cohomology wavefunctions, they do not have to be the 
same in general. For example, the Levin-Gu wavefunc¬ 
tion does not include this complex conjugation, and in 
that case T 2 = T* ^ tI. In that case, the bulk and edge 
theories are therefore different: as discussed below, while 
the edge theory (TiTj^) describes a free boson with c = 1 , 
the bulk theory (TiT^) is a highly non-trivial c = — 7 non¬ 
unitary, logarithmic CFT. A similar situation was shown 
to arise for the Haldane-Rezayi Fractional Quantum Hall 
state HO]. 

It is interesting to note that the transfer matrix of the 
edge theory is always Hermitian, as it should be if it is 
to describe the edge dynamics, while the transfer ma¬ 
trix for the bulk theory can be non-Hermitian (and the 
corresponding CFT non-unitary). Note that in this sec¬ 
tion, “edge theory” refers to the theory giving the entan¬ 
glement spectrum, as this is a well-defined theory for a 
given wavefunction, unlike the theory describing a phys¬ 
ical edge, which depends on the arbitrary choice of ad¬ 
missible terms added to the edge Hamiltonian. 


We now decompose the product of group cocycles into 
three factors 

’4^[(p-,<P+\= A-[(j)-]T,i>^^<i,^A+[(j)+\ (5) 

where 4)±(?/) = </>±(0*,j/) and where T is either Ti or T 2 
depending on where the cut lies (we will take T = Tj in 
the following without loss of generality and in agreement 
with Figure [^. A±[( 5 ()±] is the product of cocycles over 
all triangles strictly inside {r^}. 

Using the fact that |A+[(/)+] p = 1, the cocycle factors 
disappear for all triangles strictly inside {c’'"}, and one 
can trivially compute the sum for all sites with a; > 0 '*', 


IV. IDENTIFICATION OF CFTS 

After having established the usefulness of Z for stuy- 
ding both the bulk wavefunction and the entanglement 
spectrum, we now establish its criticality and identify the 
emerging CFT for a large set of examples. Before delving 
into a more specific analysis, several guiding principles 
should be identified. 

The first one concerns symmetries. The microscopic 
wave functions coming from group cohomology, and con¬ 
sequently also the partition function Z, all obey the sym¬ 
metry group of the SPT G when placed on closed bound¬ 
ary conditions. However there is often a large freedom 
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of choice in writing down such wave functions which can 
lead to extra on-site symmetries. As a starting point, 
and in order to facilitate the calculations, it is natural 
to make the most symmetric choice available. In prac¬ 
tice we find that this amounts to choosing a particular 
branching structure with an hexagonal unit cell. 

The second one concerns the relation to a classical Sta¬ 
tistical Mechanics (Stat. Mech.) model. The partitions 
functions Z obtained here all have complex Boltzmann 
weights. However, in all examples given below (except 
the Levin-Gu wavefunction and the case of G = Ga with 
p 3), we find that judiciously summing-out a particular 
sublattice of the triangular lattice leads us to a true sta¬ 
tistical mechanical models with real, positive Boltzmann 
weights. These classical models appear either as loop 
models or spin-ice like models having zero discrete di¬ 
vergence constraints. Besides being amenable to Monte 
Carlo simulations, in many cases these models can be 
solved by promoting the discrete degrees of freedom to 
continuous ones and establishing a suitable map to a com¬ 
pact boson. 

Remarkably we find three classical models obtained 
from group-cohomology wave functions which are in fact 
integrable models for which the CFT can be obtained in 
a rigorous manner. In several other cases where we are 
not aware of mappings to integrable models, our numeri¬ 
cal results show some hints of integrability: CFT-implied 
degeneracies are exact on the lattice and finite size cor¬ 
rections appear particularly small. This raises the in¬ 
triguing possibility that group cohomology cocycles are 
linked with integrable models. 


A. Abelian Symmetries 

1. Type i cocycles 

First is the group cohomology wavefunctions for an 
SPT with a. Zm symmetry [2]. To write such wavefunc¬ 
tions, one must choose a branching structure, or equiva¬ 
lently, a consistent ordering of the vertices, (ri, r 2 , r^), of 
each triangle on the lattice. Crucially for later analysis, 
we choose this ordering such that, for each triangle, the 
first, second and third vertex belongs to the A, B and C 
sublattice, respectively (see Fig. [^. In the following, we 
will refer to this choice as the ABC branching structure. 
In 2D, there are N different paramagnetic phases with 
Zjq symmetry. If we index them by p = 0,..., A — 1, a 
wavefunction belonging to the p-th SPT phase is given 
by 


(ri,r2,r3)GT 

f{4>l 4>2 4>3) = ^*^2 - <(>l) + (03 - ^ 2 ) - {(1)3 - <('l) 

^ (7) 



FIG. 2: Advantageous branching structure for which points 
on the A sublattice (shown in blue) always appear first in 
the ordering of vertices in their surrounding triangles. The 
B and C sublattices are shown in red. In each triangle, the 
first vertex has two outgoing edges, the second vertex has one 
outgoing and one incoming edge, and the third vertex has two 
incoming edges. This branching structure is referred to as the 
ABC branching structure in the following. 


where T is the set of triangles with ordered vertices in¬ 
dexed by (ri, r 2 , r^), S(ri,r 2 ,r 3 ) is 1 (~1) if tbe vertices are 
clockwise (anticlockwise) oriented, i/p is known as a cocy¬ 
cle, and (^ —</>'} = (((/> — (/)'+ 'iN/2\ ) mod N) — [A/2J. 
This choice of cocycle is a slight modification of the more 
common choice sn, where no addition or substraction of 
[A/2J is required. 

A technical advantage of this branching structure is 
that any A sublattice point, say ra, appears first in the 
ordering of the vertices in all the triangles surrounding 
it. Consequently the ratio of over A^^^y, where 

{(j)r} + lr„ is the set {(()(,} with = cj)^ + Sr^r^ mod N, 
depends only on the points surrounding (indexed in a 
clockwise fashion by ra,k with k = 1,..., 6): 




b2,ra 


N 


( 8 ) 




• (9) 


For “smooth” configurations, where \ {(p — 4>') \ <^ N for 
neighboring sites, the function d 2 ,r„ G Z is a discrete 
analog of the notion of vorticity [25]. Since d 2 ,r„ =0 
on the torus, we obtained an extra global Zjv symmetry, 
corresponding to rotating (p only on the A sublattice. 

The partition function Z can now be simplified con¬ 
siderably by summing out all the variables on the A sub¬ 
lattice {(pr^)- For N and p coprime, the cancellation of 
phases of the exponential in Eq. leaves only the sub¬ 
set of configurations on the B and G sublattices {prb,a}o 
for which d 2 ,ra mod A = 0 on every A sublattice point. 
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This yields 

z = E E = "*-* E 1 (10) 

W'-t.clWr-al 

where is the number of A sublattice sites and we 
also used the cocycle property that r'm(0,0, (/) 2 , = 

1 V (j)2,4’3- Notably, we arrived at a classical partition 
function with real positive Boltzmann weights and thus 
any CFT that might describe local observables in this 
theory should be unitary. 

To obtain the CFT associated with Z, we follow 
a heuristic approach which works for spin-ice models 
wherein a discrete, lattice zero-divergence constraint is 
promoted to an equation of motion for a continuous held 
[32] ■ The analog of the spins in the spin-ice model is given 
by the link variables {(j)i — (j)j) of our model. A slight dif¬ 
ference is that d 2 (ro) = 0 is a zero-curl constraint for 
these link variables rather than a zero-divergence con¬ 
straint. Consider promoting ^ to a compact boson ip 
using the embedding p = 2 'e<^IN G [0,27r), a possi¬ 
ble candidate for the Lagrangian is £ = Us¬ 

ing a duality transformation, the Lagrangian becomes 
£ = 4 ^(V 0 )^, with Q being the dual held of p> such 
that is the vortex density of p m- The equation 
of motion of 6 is given by = 0 and can be seen as 
a continuum version of the zero-curl constraint we have 
found on the lattice. 

Based on the above approach, a sensible continuum 
action is given by 

^ = / + m=N + 

+ y{Vm=N + Vrn=-N) 

where g is the stiffness, Ue.m is the vertex operator of 
electric charge e and magnetic charge m, with scaling 
dimension /S.e,m = ^1^9 + 5 w^/ 2 , and where only the 
most relevant vertex terms allowed by symmetry were 
kept. In terms of p and its dual variable 0, the elec¬ 
tric and magnetic operators are given by 14 = and 
Vm = e*"*®. Considering “smooth” conhgurations of </>, 
the microscopic d 2 ,r„ = 0 condition translates under the 
embedding to allowing only magnetic charges which are 
multiples of N, i.e. m G NZ, in the action. Likewise, the 
discrete nature of (p restricts the electric charges allowed 
in the action to respect the subgroup of U{1), i.e. 
e G NZ. If both these magnetic and electric terms are 
non-relevant, which is true for A^^/4 > g > 4/iV^, the 
theory is attracted to a Gaussian fixed point with c = I 

m- 

To associate CFT primary operators with the 0±{r) 
appearing in Eq. [^ we require two reasonable assump¬ 
tions in addition to the above embedding: (1) The CFT 
operators should transform according to a given represen¬ 
tation of the sublattice symmetry group. In the present 
case, after having integrated out the A sublattice, the 
only remaining sublattice symmetry corresponds to an 


interchange of the B and C sublattices. (2) The CFT 
fusion rules should be consistent with the Z]\j additive re¬ 
lation for the charges. This is made possible by the pres¬ 
ence of the vertex terms in the action, which corresponds 
to a finite density of screening charges with e, to G iVZ. 

For 0± (r) operators on the B or C sublattice, the em¬ 
bedding implies ^ ^ yielding 0± -G 

Ve=^i. Numerically (see Appendix [A|, we find that there 
is a sign difference between B and C: on the B (C) sub¬ 
lattice, we have 0± -G Ve==Fi (0± —14=:fi), implying 

that Ve=zfi picks up a minus sign under the sublattice 
symmetry, as allowed by assumption (1). For a charge 
a on the A sublattice point r^, the summation over 
enforces a non-zero vorticity d 2 ,ra whose value is deter¬ 
mined by pd 2 ,r„ — a = 0 mod N. The embedding thus 
associates a non-zero discrete vorticity d 2 ,r„ for ^ to a 
magnetic charge with to = d 2 ,r„ for p. For p = I, this 
simply leads to 0± —>■ I4i=±i- For p > I, the discus¬ 
sion remains the same if, on the A sublattice, one enu¬ 
merates the charge basis by the position of charge ±p 
particles instead of charge ±1 particles, since one has 

—>■ I4t=±i. Different values of p therefore simply cor¬ 
respond to a reshuffling of the different charge values, and 
we will use p = I in the following in order to simplify no¬ 
tations. In short, charges on the A sublattice correspond 
to magnetic charges and charges on the B and C sub¬ 
lattices correspond to electric charges (see Fig. |^. The 
case N = 2, made special by the fact that 0+ = 0_, is 
discussed in Appendix |Aj 

In Appendix [^ we numerically confirmed the above 
predictions for N = 2 and = 3 in terms of central 
charge, operator content, power law behavior of correla¬ 
tion functions and representation of lattice operators. We 
find p = 1 ± 5 X 10-4 for fV = 2 and g = 0.925 ± 0.01 for 
TV = 3. Note that the value of g is in general not univer¬ 
sal, except possibly when an extra symmetry is present, 
as is the case for = 2 for which adding 1 to all the pr 
on any given sublattice is a symmetry. This symmetry 
between the three sublattices corresponds to a symmetry 
between electric and magnetic charges in the CFT lan¬ 
guage and thereby forces g to be at its electric-magnetic 
dual value of 1 |4^ . 

It is now illuminating to study the symmetry action 
at the light of the previous results. The global symme¬ 
try action corresponds to applying p ^ p 1 on each 
site. At the CFT level, this symmetry action translates 
into different transformations depending on the sublat¬ 
tice. On B and C, the symmetry action translates via 
the embedding to p —>■ p -I- On the A sublattice, in- 

. ■■ . . —27ript^ 

serting the microscopic operator e « corresponds to 
adding a vorticity of one, or equivalently to inserting a 
Vm=i = e*® operator in the CFT language. It is then nat¬ 
ural to associate the microscopic variable —2'Kpp/N on 
the A sublattice with 9 in the CFT and, consequently, to 
associate the microscopic symmetry action p ^ p +1 to 
9^9— Thus, the global Zisi ^2 symmetry, for non¬ 
zero p, rotates both 9 and p simultaneously. The above 
symmetry action is exactly the one discussed in Ref. |TT| 
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with regards to the edge theory of such SPTs. As shown 
in that work, given such an action of the symmetry, the 
only way a symmetry-preserving term can gap out the 
theory is by inducing a spontaneous symmetry breaking. 

We now go back to the microscopic derivation of 
the wavefunction amplitude in the symmetry-charge 
basis, A{±. j..}. Focusing on > 2, we remove 

fast sublattice oscillations by defining = 

where Sr-^c is 1 (0) if is on the 


C {A, B) sublattice, and the set of electric and mag¬ 
netic charges {{ei,mi,ri}) is determined from {±,ri} 
by the previous association (e.g. (+,rQ) —> ( 0 ,-fl, ra)). 
The wavefunction amplitude can now be written as 
ei,mi{ri)), where the average value 


is taken with the action given by Eq. 11 The usual 


renormalization group procedure can be used from the 
lattice scale Og to the scale of average interparticle dis¬ 
tance I [p = leading to a suppression of A (A') as 

^ with A = N"^I2g (A = gN‘^/2). We focus 

first on the infinitely dilute charge limit p —> 0 , for which 
A, X' —>■ 0, yielding a simple form for the wavefunction 
amplitude: 



FIG. 3: Electric particles live on the B and C sublattices 
(shown in red) and magnetic particles live on the A sublattice 
(shown in blue). The blue curvy arrows indicate the magnetic- 
particles-generated flux that is seen by the electric particles. 


A 


{ei,mi,ri} 


k<lGE s<tGM 


n 


/ zk - zt 

Vkfe - zt\) 


( 12 ) 

where z = +i Ty and where we used the fact that each 
particle is either electric or magnetic, but not both, to 
separate the particle indices into two subsets, £ and Al, 
for electric and magnetic charges, respectively. Notably, 
unless the above correlator vanishes. 

Consequently, we find an enhanced C/(l) x 17(1) symme¬ 
try compared to the original wavefunction for which the 
above neutrality condition was only obeyed modulo N. 

At small but hnite p, A and A' are non-zero and act 
as a condensate of charge ±N particles thereby breaking 
the U{1) X 17(1) symmetry present at p — 0 down to the 
initial Zjv x symmetry. This mechanism is explained 
in Appendix where we also give a self consistent ar¬ 
gument showing that the Laughlin picture persists even 
with a finite condensate density. 

Interestingly, in Ref. [45] a similar 17(1) x 17(1) 
wavefunction was proposed (however with only posi¬ 
tive charges). Consistently with our picture, condensing 
charge N particles at the level of the Chern Simons the¬ 
ory associated with that wavefunction was shown to yield 
the N different Zj^ SPT phases [llj . 


2. Type ii cocycles 

Let $ = ( 01 , 0ii) denote an element of Z^ x Z^, with 
N prime. In the following, we will use p = I, II as the 
index that distinguishes between the two groups. Accord¬ 
ing to Ref. (SSI), the third cohomology group oi Z^ x Zjsi 


TABLE I: Association of GET primaries and microscopic op¬ 
erators for several diffent cocycles of G = ZnX Zn ■ Regarding 
microscopic operators, 0 a ,77 stands for the insertion of an op- 
erator e* « on the A sublattice and 4>b(c),r] stands for the 
insertion of an operator e'^r^-n on the B or C sublattice. The 
GFT operator Ve(m),ri is an electric (magnetic) operator for 
the compactified boson indexed by g. 


Pi 

Pii 

p' 

V^,i 

Ve,! 

Vm,U 

Fe,II 

Gentral charge 

AO 

0 

0 

0 a,I 

0 ij(c),I 



1 

AO 

AO 

0 

0 a,I 

0 i>(c),I 

0 a,II 

0 i>(c),II 

2 

0 

0 

AO 



0 a,I 

0 i 7 (c),II 

1 

0 

AO 

AO 



0 a,I + 0 a,II 

0 i>(c),II 

1 


is spanned by products of two distinct types of 3-cocycles. 
The type i cocycles involve just one of the Z^ groups 
while type ii cocycles involve both. Namely, 


, _ _ _ ^ 27ripTj01,77/77 

2-Kip' (pi ifii 


j, _ ( 02,77 ~ 4'l,v) + {4>3,rj — 02 , 7 /) — (03,7/ — 01 , 77 ) 

h- ^ 

( 13 ) 

where we again choose the branching structure given in 
Fig. such that, for each triangle, <i)i (resp. 4 ) 2 , $ 3 ) 
refers to the site on the A (resp. B, C) sublattice. A 
generic cocycle is then labelled by (pi,pii,p') and given 

l)y ■ 

We consider first some simple cases, as summarised in 
TablejH First, when the only non-trivial cocycle is a type 
i cocycle for one of the two groups, say the first one, the 
previous arguments for G = Z^ can be applied to the 

































7 


first group, while the degrees of freedom of the second 
group form a trivial paramagnet. Second, when there 
are two non-trivial type i cocycle, one for each group, one 
finds two decoupled G = Zm SPT phases, and again the 
previous arguments apply directly. Third, when p' ^ 0 
but Pi = pii = 0, one should sum (pi on the A sublattice 
and obtain a compactihed boson associated with the 
coarse graining of pu degrees of freedom on the B and 
C sublattices. For N > 2, the electric charge operators 
Ve=±i of that boson would correspond to the microscopic 
operators e « on the B and C sublattice as before. 
However, the magnetic charge operators Vm=±i would 

now correspond to e « on the A sublattice. All the 
other microscopic degrees of freedom {pi on the B and 
C sublattice and pu on the A sublattice) form a trivial 
paramagnet. Lastly, for p' ^ 0 and pu ^ 0 but pi = 0, 
one obtains a single compactihed boson ipn associated 

^ i 27r2(i) j j 

with pu. For N > 2, the microscopic operators e w 
on B and C correspond to 14=^1 and the microscoic op¬ 
erator e ^ on the A sublattice corresponds to 

Vm=±l- 


A more complicated case, which does not follow di¬ 
rectly from the previous G = results, \s p' ^ Q 
and Pi ^ 0 but pu = 0. In this case, one can in¬ 
tegrate out pi on the A sublattice to obtain that the 
condition (pid 2 ,r„[ 4 i] A mod iV = 0 has 

to be satished for each B,G hexagon centered at ra- 
The discrete vorticity d 2 ,r„[ 4 ] is dehned in Eq. It 
is then natural to consider the linear transformation 
{Pi,Pu) = {PiPhPiPu) + {p'Pii,-p'Pi), dehned in the 
two-dimensional vector space over the hnite held Fjv 
(which is well-dehned since N is prime). This linear 
transformation is an isomorphism if its determinant is 
non-zero, i.e. if (pj +p''^) mod N 0. We now focus on 
the case A^ ^ 1 and pi,p' A I for which this condition is 
fulhlled. 


One can then repeat the assumption made previ¬ 
ously that the important conhgurations are those in 
which pi and pu change by much less than N/2 be¬ 
tween adjacent lattice sites. In this case, one has that 
Pid 2 ,ra [Pi\ + yd 2 ,r-„ [4ii] = ‘^ 2 ,ra [4i] ^md a compactihed 
boson pi emerges, corresponding to the coarse graining 
of pi- Repeating the rationale used for G = Z^, one 

^ , , . . ■^ 2 Tri 4 >j 

would then hnd that the microscopic operator e « 
on the A sublattice corresponds to I4„=±i and that the 

±27r2<^I 

microscopic operator e « on the B or C sublattice 
corresponds to I4=±i- All the other microscopic degrees 
of freedom {pu on A and pu on B and G) form a trivial 
paramagnet. A test of this conjecture, as well as a gen¬ 
eralization for all (pi,Pii,p 4 S'Hd non-prime N, is left for 
future work. 


3. Type 3 cocycle and exact mapping to loop models 

When a product of more than two Zn groups is consid¬ 
ered, the third cohomology group is spanned by the type 
i and ii 3-cocycles considered above, as well as one ad¬ 
ditional distinct 3-cocycle, called type iii [15]. Focusing 
on the case oi G = Zn x Zjq x Z^ and denoting group 
elements as $ = {pi, pu, Pin), the type iii 3-cocycle is 
given by 


/n ^ ^ X. \ 27ripaibnein 

i"iii,p( 0 ,«'i,«' 2 ,$ 3 ) = e « (14) 

a-i = Pi,i 

bu = ( 42.11 — 4 l,Il) 

Cm = (43,111 — 42,III) 


where we again choose the branching structure given in 
Fig. such that, for each triangle, ihi (resp. 4 ) 2 , $ 3 ) 
refers to the site on the A (resp. B, C) sublattice. No¬ 
tably this cocycle is not defined in terms of the / factor 
present in the previous examples (Eqs. and 13) and 
the reasoning used for G = Zjq appears irrelevant in this 


case. 


Instead, consider re-writing this 3-cocycle as 


t'iii,p(0, 4>i, 4>2, 4>3) = i7p(0, $ 1 , $ 2 ) 


(15) 


where t&j = (4i,i) 4i,ii) and where Vp can be recognized as 
a 2-cocycle associated with aG = Z^x Z^ symmetry [53] 
(i.e. a cocycle in the second cohomology group ^{^{Zn x 
Zn,U{1))). 

This suggests that type iii 3-cocycles are described by 
the decorated domain wall picture of SPTs gTj. Follow¬ 
ing this picture, each domain wall for the 4iii degrees 
of freedom should be thought of as carrying its own G = 
Zjq X Zjq ID-SPT phase for the 4i and 4ii degrees of free¬ 
dom. [As a technical comment, note that the 2-cocycles 
of the form given by Vp span BA {Z^ xZn,'H^{Zn, U{!))), 
which is the subgroup of TL^{Zi^ x Z^ x Zf.j,U{l)) for 
which the decorated domain wall picture applies mi 

The N = 2 case: For N = 2, the only non-trivial case 
corresponds to p = 1. Furthermore, since the cocycle is 
real in this case, one has lymu = and one can ignore 
the orientation of the triangles. First, let us consider 
the restriction of 4iii {r) to its values on the C sublattice 
only, 4iii(^c)- Since the G sublattice is itself a triangular 
lattice, its dual lattice is an hexagonal lattice, and it is 
the one formed by the A and B sublattices. As shown 
in Fig. 1^ this means the domain walls of 4iii(^c) form 
non-intersecting loops on the A, B hexagonal lattice. We 
now show that these loops carry ID SPTs for (4i, 4ii)- 

From Eq. |T^ it is clear that the only triangles with a 
non-trivial factor are the ones for which cm = 1, i.e. the 
ones for which 4iii is flipped when going from the B site 
to the G site of the triangle. If this is the case, we call the 
triangle active and its non-trivial factor is given by the 
2 -cocycle Vp for the (4i, 4ii) degrees of freedom living on 
the two sites of the AB edge of the triangle. Now, since 













FIG. 4: Decorated domain walls picture emerging naturally 
for our choice of cocycles and branching structure for the case 
oi G = Z 2 x Z 2 X Z 2 - The arrows give the value of (f)iu{rc) 
(the restriction of (j>m to the C sublattice) and the shaded 
blue lines give the domain walls of the 4’ii\{rc) configuration, 
which form non-intersecting loops on the hexagonal lattice 
formed by the A and B sites. AB edges are shown in blue, 
BC edges in red and AC edges in green. 


each AB edge (shown in blue in Fig. is shared by two 
triangles, there are two possibilities: either (1) the AB 
edge is occupied by a (l>iii{rc) domain wall (highlighted 
in shaded blue in Fig. |^, in which case the C sites of 
the two triangles have a different value for (pm, or (2) 
the AB edge is unoccupied by a (puiifc) domain wall, 
in which case the C sites of the two triangles have the 
same value for pm- Crucially, for the first possibility, 
we know that one of the two triangles is active and the 
other one is inactive, since, when going from one C site 
to the other along the BC edges of these two triangles, 
there has to be exactly one flip. In this case, one 
obtains a factor of Vp for the AB edge. For the second 
possibility, however, we know that the two triangles are 
both inactive or both active, since, when going from one 
C site to the other along the BC edges of these triangles, 
there has to be, respectively, zero or two pm flips. In this 
case, the overall factor is trivial, since = 1. 

Now that we have shown that a Vp factor appears for 
each AB edge covered by a pniirc) domain wall, we turn 
our attention to the {pi,pii) degrees of freedom along 
one such domain wall. Since the vertices along this do¬ 
main wall alternate between the A and B sublattice, it is 
advantageous to split the domain wall as a succession of 
AB pairs indexed by i and for which the two vertices in¬ 
side the pair are indexed hy u = A, B. This leads to the 
following notation for the degrees of freedom along one 
domain wall: {pi^u,iTpi,u,ii)- The Izp factors then generate 
the following overall phase factor for the entire domain 
wall: 


— g®'"' —0i-l,B,Il)] 

= g“ Z]i[0i,B,Il(0i,A,I-0i + l,.4,l)] 


(16) 


where we have used the Z 2 algebra obeyed by p and 
the last line follows from discrete integration by parts 
(and assuming a closed domain wall, which is always the 
case for periodic boundary conditions). As shown else¬ 
where na Ea, these phase factors describe a ID SPT 
ground state with G = Z 2 x Z 2 . Furthermore, correla¬ 
tions in the ID Stat. Mech. model defined by these am¬ 
plitudes are long ranged [30] ■ Indeed, one may integrate 
out all pi^A,i (resp. pi,B,ii), and obtain that all pi,B,ii 
(resp. pi,A,i) must be equal, and hence are long-range 
correlated. 

Collecting these results, we find that there is a two to 
one mapping between pm configurations and loops (i.e. 
domain walls) on the hexagonal AB lattice. For a given 
pm configuration, the correlation for the remaining de¬ 
grees of freedom {pi, pu) at different sites is strictly zero 
unless these sites lie on the same pm{rc) domain wall, 
in which case the correlation is infinitely-ranged along 
the domain wall. Calculating a correlation of the type 
Gii(rb,r(,) = therefore amounts to eval¬ 

uating the probability of rj, and r(, lying on the same 
Piiiifc) domain wall. This is precisely the definition of 
the correlation function of 2 -leg watermelon operators, 
as defined in Section UV Cl 

To determine the behavior of the above watermelon op¬ 
erator, we need to identify what is the precise loop model 
that was obtained. These are defined by two parameters, 
n and x, denoting the loop fugacity and loop tension, 
respectively. Fixing pm, and summing up, say, pi, one 
obtains that pu degrees of freedom along a loop are all 
identical. Since they can all either be 0 S ^’2 or I G ^ 2 , 
this gives a weight of n = 2 for each loop. Loop ten¬ 
sion comes from the fact that the degrees of freedom 
that lie on loops cannot fluctuate freely. Thus the pres¬ 
ence of a loop of length I reduces the partition function 
by a factor of 2“*/^. The loop tension is therefore given 
by a: = I/-\/2- Interestingly, this loop model is right at 
an exactly solvable point of the honeycomb seven-vertex 
model |48| described by the SU{2)i Wess-Zumino-Witten 
(WZW) CFT. In terms of the Coulomb gas description 
given in Section |IV C[ this corresponds to having zero 
background charge (cq = 0 ), and this theory is therefore 
a free boson with 5 = 1 and c = 1. The probability that 
two points at a distance r are on the same domain wall 
is given by = r~^, where A; = — 60/25 -|- gP/2 

gives the scaling dimension of the 2 Z-leg watermelon op¬ 
erator. This exponent governs correlations functions of 
microscopic charge operators of the type 6 *’^'^''“/'’’^^", as 
exemplified by Cii{rb,r{). 

The N = 3 case: Considering the case of iV > 2, 
domains walls may now branch and intersect. Repeating 
the previous reasoning, one ends up with loop models 
in which the loops may branch. More accurately, each 
domain-wall carries a group element and the branchings 
must obey the loop algebra. Such loop models cannot 
be solved straightforwardly using the Coulomb gas ap¬ 
proach and thus we resort to a numerical calculation of 
the transfer matrix for N = 3. We found again a critical 
















9 


model, with a central charge of c = 2±0.05. The first ex¬ 
cited state of the transfer matrix corresponds to a scaling 
dimension of 0.4 ± 0.008 and is 6 -fold degenerate. Given 
the abundance of compact bosons we have found so far, 
a natural guess is that the CFT obtained here is that of 
two compact bosons. However, if these are decoupled, 
it is impossible to reproduce the observed degeneracy of 
6 for the lowest excited state using just the two param¬ 
eters gi,g 2 - Considering two coupled compact bosons 
[26] . the moduli-space (or the space of critical theories) 
is spanned by four parameters [55] ([Gn, G 12 , G 22 , . 612 ])- 
This space includes the SU{3)i WZW model model at 
[1,0.5,1, 0.5] and so all these models can be thought of as 
exactly marginal deformation of this model. Note that we 
have already encountered similar deformations since the 
models obtained previously for G = can be viewed 
as exactly marginal deformations of an SU{2)i WZW 
model. While the central charge and 6 -fold degeneracy 
can be reproduced by choosing the four aforementioned 
parameters within the appropriate region of parameter 
space, our finite-size numerics do not allow us to make a 
definite statement about the identification of this theory 
with a two-component Luttinger liquid. 

To sum up, the G = [Zn= 2 ]^ type 3 SPT corre¬ 
sponds to an SU{ 2)1 WZW theory with c = 1, while 
the G = type 3 SPT corresponds to a critical 

theory with c = 2 that could be related to a marginal de¬ 
formation of an SU{3)i WZW theory. This qualitative 
difference between N = 2 and iV = 3 is highly unusual 
and seems to indicate that the usual Chern-Simons rea¬ 
soning m within which one embeds the finite group in 
a product of U{1) group does not apply in a straightfor¬ 
ward way to type 3 cocycles. Interestingly, the non-linear 
sigma model construction m embeds the finite group 
in SU{2) X SU{2) for I? = 2 spatial dimensions. While 
[Zjsr= 2 Y is afinite group of SU{2)xSU{2), [Zjq^^Y i® i^ot, 
and a qualitative difference between these two cases is 
therefore consistent with this picture. It would be highly 
interesting to study what happens for N > 3. 


B. Non abelian symmetries: D 3 

The simplest non-Abelian group to consider is the 
Dihedral group of the triangle (D 3 ). It has six el¬ 
ements {g = (4>, 0 )) corresponding to three rotations 
{9 G [0,1,2]) and an inversion $ G [0, Ij. Due to the non- 
Abelian nature of the group, it makes sense to use multi¬ 
plicative notation for the group operation. Multiplication 
is given by = (($ ± $') 2 , ((- 1 )^'^ ± 

where ■ ■)j~^ stands for modulo N. The cohomology 
group here is 'H?{D 3 ,U( 1 )) = Z^ and so 6 distinct SPT 
phases exist (including the trivial phase). The cocycles 
can be written as (see 6.20 in Ref. 021) 

a;p(($i,0i),(ci>2,02),(4>3,03))= (17) 

^{_^J(_^^*2 + f3ei[(_l)t'3e2+e3-((-l)'*'3e2+e3)3]+|$i$2$3]} 
= Vp (1, ($l,0l), (4 >i,6»i)($2,6'2), ($i,6'i)(4>2,6»2)($3,03)) 


For p = 3, a certain simplification arises: the factor 
[(_ 1 ) 3 > 302-|-03 — ((— 1)^302 ±^* 3 ) 3 ] is always a multiple of 
3 and so it is effectively zero inside the above exponent. 
Consequently, the cocycle becomes which is 

no other than the type 3 cocycle for G' = Z 2 x Z 2 x Z 2 
defined in Eq. after the following identification: on 
the A (resp B, C) sublattice, $ maps to (j)i (resp (/)n, 
(/)iii). Thus, up to a multiplication by some trivial factor 
coming from the completely uncorrelated 9 degrees of 
freedom, the partition function for p = 3 and for the type 
iii cocycle are identical. They are thus both described by 
a c = 1 , g = 1 compact boson. 

For p = 1, we turn again to numerics. The large size 
of the group reduces our ability to approach the scaling 
limit. Still, looking at a circumference of L = 6 , we 
find that the leading scaling dimensions agree very well 
with c = 1, 5 = 1 (see Table (IV)). This suggests that 
for all values of p, we obtain a c = 1 , p = 1 compact 
boson. Notably, however, the scaling dimensions here 
do not show exact degeneracies on the lattice, as was 
the case for all the Abelian cases. We attribute this to 
the fact that, for p = 1, even with the ABC branching 
structure, the wavefunction does not have an enhanced 
D 3 X D 3 symmetry. 


C. The curious case of the Z 2 Levin-Gu 
wavefunction 

A different type of wave function to consider is the 
Levin-Gu wavefunction |3] for the Ising (G = Z 2 ) SPT 
whose phase factors are not directly given by group- 
cohomology cocycles. For clarity, we consider (j)^ = 0,1 
as denoting the = ±1,-1 eigenvalues of a spin 
1/2. Correspondingly, the charge = 0,1 denotes the 
= ±1,-1 eigenvalues. The amplitude A^^zy is sim¬ 
ply given by (—1) to the number of domain walls in the 
cr/ configuration. This model flows to a stable critical 
phase with central charge c = —7 described in terms of 
a Coulomb gas with background charge eo = 2/3 and 
stiffness g = 1 — Cq [301 SHI [5DH53] . There are two 
types of primary operators in this theory: electric op¬ 
erators Ve{r) = with e G Z and with scaling di¬ 

mension Ae = e(e — 2 eo)/ 2 g and 2Z-leg “watermelon” 
operators H/(r) with I G Z and with scaling dimension 
A; = —eQl2g + gl^l2. The latter correspond to imposing 
21 legs of domain walls emanating or closing at a given 
point [48| . 

Numerically, we find a good agreement with Coulomb 
gas predictions, despite the fact that we are dealing with 
a logarithmic CFT |S3]. As explained in Appendix]^ the 
same prescription as before is used to find a correspon¬ 
dence between CFT primary operators and the micro¬ 
scopic operator 0(r) — In this case, we find 

that these operators lead to some linear combination of 
14=1 and of the 2-leg watermelon operator. Interestingly, 
the 2-leg watermelon operator has the same scaling di¬ 
mension as 14=1 and the 4-leg watermelon operator has 
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the same scaling dimension as 14 = 0 , consistently with a 
Z 2 charge interpretation. We comment that, following a 
suitable choice of basis, the wavefunction amplitude still 
appears as a single CFT correlator (see Appendix . 

The attentive readers will have noticed that this is the 
second non-trivial G = Zi SPT wavefunction we have 
discussed in this work. The first one was the group 
cohomology wavefunction, and its bulk CFT was given 
by a free boson (c = 1). It is interesting to note that 
these two wavefunctions have the same entanglement 
spectrum [78] but different bulk CFTs. While we leave 
for future work the discussion of the possible implica¬ 
tions of this different bulk CFT, we note that the local 
unitary transformation taking the Levin-Gu wavefunc¬ 
tion to the group cohomology wavefunction can be writ¬ 
ten as Hr < ^ and was shown to have a 

non-trivial action on weak SPT indices [54| . Another 
perhaps relevant observation is the fact that the central 
charges of the two different CFTs differ by 8 , which could 
be an indication that the two states differ by a bosonic 
Es X Es bosonic SPT (with central charge 8 and chiral 
central charge 0). The physical Z 2 symmetry would then 
have to act non-trivially on the Eg, x Eg state in order to 
protect it. 


Considering the correlator C'(r) = 

we obtain 

_ |r|““ ^2 

{±,.n} ( 21 ) 

(X |r|““ for |r| —>• 00 

where 0 and r are on the A {B or C) sublattice, where 
a = gl2 {a = l/2g) and where .F[{±i, rj; (-I-, 0), (-I-, r)] 
is the energy of a plasma as defined in Eq. but with 
charges of magnitude -s/S and with two additional test 
particles of magnetic charge ±l /-\/2 (of electric charge 
=F1/-\/2, respectively) at positions 0 and r. In order to 
obtain the second line of Eq. one simply has to assume 
that the magnetic (electric) plasma screens. The screen¬ 
ing condition for such plasmas is g < 2 {g > 1/2) [57ll61j . 
Thus, using flux attachment, we find that |^) possesses 
quasi-long-range hidden order associated with electric op¬ 
erators which break the Z^ symmetry. 

The Laughlin-like picture put forward here can also be 
used to study the charge fractionalization appearing after 
the insertion of symmetry fluxes (see App. E- 

VI. INTEGRABILITY 


V. HIDDEN ORDER 

As we showed in the previous section, the Z^ x Z^ 
wavefunction can be written as 


A 


{ei,mi,ri} 




Zs - Zt\ 


k<ieS 


n 

k£e,t£M 


Zk - Zt 
\zk - Zt\ 


s<tGM 

ekmt 


(18) 

Using this formula, we can discuss the presence of 
hidden order. Based on this Laughlin-like picture, one 
can readily unveil a composite-particle [55] or hidden or¬ 
der [5^ structure behind these SPTs. Following Ref. jSS], 
we use a suitable flux attachment transformation (U) to 
remove complex phases from |'(/)), so that the following 
wavefunction is obtained 


= E (19) 




From Eq. (18) we find that the resulting amplitudes 


|A{ei,mi.ri}I are related to two decoupled (electric and 
magnetic) two-component plasmas via | = 

with 


E[{ei,m*,rJ] = - E — log(|rfe - nl) 


- E WsTOtglog(|rs - rtl). 


( 20 ) 


s<teM 


In this section, we discuss the relation between SPT 
wavefunctions coming from the group cohomology con¬ 
struction and integrable models. We show that (1) in 
several cases, there is an exact mapping between the aux¬ 
iliary Stat. Mech. model of these SPTs and an integrable 
model, ( 2 ) even when such a mapping could not be found, 
the numerically obtained finite-size spectrum exhibits the 
exact degeneracies predicted by CFT and (3) the finite- 
size corrections to scaling dimensions show better con¬ 
vergence than what would generically be expected. 

We have found three examples where the auxiliary 
Stat. Mech. models corresponding to a group coho¬ 
mology SPT wavefunction can be exactly mapped to 
an integrable model. The first example is the case of 
G = Z 2 X Z 2 X Z 2 with the type 3 cocycle, for which we 
have established an exact mapping to a loop model with 
loop fugacity n = 2 and loop tension x = 1/ v^, which 
happens to be an integrable point. At these values, it is 
known that its scaling limit is given hy a c = 1, g = 1 
compact boson, or in other words a SU{2)i WZW the¬ 
ory [35]. The second example appears for G = Z 2 and 
is obtained in a similar way since the partition function 
of the case oi G = Z 2 x Z 2 x Z 2 is equal to the one of 
G = Z 2 times that of a trivial paramagnet. The mapping 
between the two descriptions works in the following way: 
on the A (resp B, G) sublattice, (j) maps to <j)i (resp (pu, 
t/ni)- The third example is that of Dg discussed in Sec. 
11X11 which again maps to the G = Z 2 x Z 2 x Z 2 case. 

Considering the case of G = Z 3 , it cannot be mapped 
to the same n = 2 loop model since the numerically ob¬ 
tained value of g differs from 1. Still, the numerical re¬ 
sults point to some fine-tuned point: First, the numer- 
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ically obtained finite-size spectra have the exact degen¬ 
eracies predicted by CFT (see Table III I and, second, 
finite size corrections seem to be qualitatively smaller 
than generically expected. Since the leading correction to 
scaling dimensions goes like 0{1 /where A is the 
scaling dimension of the most relevant operator allowed 
in the action, a microscopic model generically exhibits 
corrections that scale like 0(\/LF‘) because of the opera¬ 
tor L_ 2 T- 2 l [Ln/Ln is the left/right Virasoro generator 
and 1 is the identity operator) [62] ■ If the lattice is trian¬ 
gular, non-conformal-scalar terms such as L~^l are also 
generically allowed in the actionjSlj, and would lead to 
0{1/L) corrections to scaling. Instead, from the numer¬ 
ics carried out for G = Z 3 , corrections to the scaling of 
the identity operator (a quantity which does not depend 
on our estimate of g) are given by 3 x 10“^ for L = 10, 
which seems to indicate that 0{1/L) corrections are fine 
tuned close to zero, and that even the usual 0{1/LF') are 
reasonably small. Similarly small corrections are present 
for all other scaling dimensions. 

This is in contrast to the case of the Levin-Gu wave- 
function, for which the loop model based on the original 
wavefunction is given byn = —1, x = 1, which itself is 
not an integrable point, but which is in the basin of at¬ 
traction of the integrable point n = — 1 , x = Ij \/2 — \/Z. 
There, numerical results obtained for the original wave- 
function show criticality but almost no similarity with the 
conformal theory predictions even for the largest system 
size we reached (L = 24). On the contrary, the numer¬ 
ics performed at the integrable point show extremeley 
good convergence to the CFT predictions, as explained 
in Sec. IIV Cl 


VII. CONCLUSION 

In conclusion, we have proposed a CFT-based method 
to analyze SPT phases which reveals their properties in 
a transparent way and we showed how this description 
emerges from a microscopic treatment n la of their 
wavefunctions. This approach allowed us to re-derive 
various expected results, such as the edge theory and 
symmetry-flux response of type ii SPTs. Concerning new 
results, it allowed us to derive a Laughlin-like picture 
for the Zj\jy 2 group-cohomology wavefunctions which in¬ 
cludes quasi-long-range hidden order and a composite- 
particle interpretation. It further provided us with the 
entanglement spectrum of all SPTs we have considered. 

From the mathematical perspective, perhaps the most 
intriguing result is the direct connection between group 
cohomology cocycles and critical, sometimes integrable, 
models. Indeed, all “Statistical Mechanics” models we 
have obtained by putting products of group cohomology 
cocycles as the Boltzman weights, were critical and some 
were even integrable. Notably, even if one associated 
these models with the physical edges, criticality is by no 
means implied as generic edge theories of SPTs can be 
both critical or broken-symmetry phases. It is possible 


that this hidden criticality of the wave function ampli¬ 
tude is generic or that group cohomology wavefunctions 
are naturally tuned to be in the critical regime where 
CFT gives us a useful handle on their properties. Re¬ 
gardless, our results suggest an intriguing correspondence 
between group cohomology cocycles and critical lattice 
models. 

It is also noteworthy that the loop model related to one 
of the SPTs we looked at has a loop fugacity of n = 2 
and therefore breaks the “n = y /2 barrier” which sets an 
upper bound on the loop fugacity of topological quantum 
loop models and has strongly hampered the realization 
of doubled S'C/(2)fc topologically ordered states with loop 
models (53] |M|. This barrier arises from the fact that 
(1) loop models are critical only for \n\ < 2[3S] and (2) 
the Stat. Mech. model considered within this construc¬ 
tion is given hy Z = where C runs over loop 

configurations. Crucially, as argued in this work, the rel¬ 
evant Stat. Mech. model for SPTs is given instead by 
Z = V’(^) therefore only has a n = 2 barrier. It 

would be interesting to see how this discussion general¬ 
ized to SPTs with G = [Zn> 2 T, in which case loops are 
allowed to branch. 

The current work suggests that a CFT approach to 
SPT phases can be unifying, useful for microscopies and 
that it can bring physical intuition about these states. 
There are various directions along which it could be fur¬ 
ther explored. It seems for example natural to gener¬ 
ate fractional (i.e. long-range entangled) SPT phases 
by orbifolding the CFT with respect to a subgroup 
G. Such orbifolding would introduce magnetic operators 
with fractional charge into the CFT which would bind, 
via the symmetry-flux response argument, a fractional 
symmetry-charge. It would also be interesting to further 
explore the mapping between SPTs and integrable mod¬ 
els introduced here, and even more so to study whether 
it may be reversed. This way, the sizable knowledge that 
was accumulated in studying integrable models could be 
used to obtain new microscopic models for SPT phases. 
Since the entanglement spectra we have obtained all had 
integer central charge, interesting candidates are inte¬ 
grable models with a fractional central charge. Further¬ 
more, since the group cohomology approach becomes less 
comprehensive for fermionic SPTs [65j . it may be useful 
to consider fermionic integrable models as a means of gen¬ 
erating microscopic wavefunctions for fermionic SPTs. 
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Appendix A: Numerical evidence for the 
identification of a free boson CFT 


1. Abelian 


Here we give numerical evidence which strongly sup¬ 
ports the conjecture made in the text, that the compact 
boson CFT is the long distance theory governing the Z 2 
and Z 3 group cohomology wavefunctions. To this end, 
we consider the statistical mechanics problem associated 
with the “Boltzmann weights”, given the hexago¬ 

nal branching structure used in the main text and the A 
sublattice analytical integration. The partition function 
(Z) then simply counts the number of configurations on 
an hexagonal lattice which obey the zero discrete vortic- 
ity condition (d 2 (ra) = 0 ). 

To analyze the CFT operator content of this problem, 
we use the transfer matrix approach |62j . We consider 
the hexagonal lattice as a brick lattice with the straight 
lines being vertical (i.e. the bricks are laid vertically). A 
transfer matrix (Tj) is then generated such that Tr[T; ' ] 
gives the partition function of a torus of circumference I 
and length L, where the torus consists of L vertical lines, 
each containing I sites. The negative logarithm of the 
eigenvalues of T; (A(^)i), when normalized according to 
A{l)i = X{l)i i^ 2 ^y( 2 y/ 3 ) shifted such that A{1)q = 0, 
then correspond directly, in the limit of large Z, to the 
scaling dimension (A^) of different operators in the CFT. 

The primary operators of a compact boson CFT are 
given by the vertex operators with e, m S Z, where 
e is the electric charge and m the magnetic charge. The 
operator content consists of these primary operators and 
of their descendants. For the microscopic model stud¬ 
ied here, we expect both electric and magnetic charges 
to be present since, as argued in the main text, charge 
operators on the A sublattice generate magnetic charges 
and charges on the B and C sublattice generate electric 
charges. With a compactification radius of 1, the scaling 
dimensions of primary operators as a function of g are 
■ Each of these vertex operators has de¬ 
scendants, or equivalently particle-hole excitations of the 
compact boson [3S]. For our analysis we shall only require 
the two lowest lying of such excitations (n = ± 1 ), with 
scaling dimension + 1- The 7 most relevant scaling 
dimensions A^, assuming g smaller than but close to 1 , 
are given by A = 0 for the identity /, A = 5/2 for the 
two unit magnetic charges Vm=± 1 , A = 1 / 2 g for the two 
unit electric charges A = 1 for the two level one 

descendants of the identity I (Si {n = ± 1 ), A = gj2 + ll2g 
for the four combinations of 14 =±i,m=±i: A = g/2 -|- 1 
for the four level one descendants of the unit magnetic 
charge Vm=±i ® (n = ±1) and lastly A = l/2g + 1 for 
the four level one descendants of the unit electric charge 


Ve=±i 0 (n ± 1). This sums up to a number of 19 op¬ 
erators for the 7 most relevant scaling dimensions. Ta¬ 
ble [TT] lists the first 19 values of A^, for the Z 2 group- 
cohomology wavefunction with the hexagonal branching 
structure, and I = 16,18, 20. The degeneracy of eigen¬ 
values is given in the brackets, when it is different from 
1. This table agrees well with g = 1, with discrepan¬ 
cies of the order of 10 “^ which are attributed to finite 
size effects. Table |III| lists the first 19 eigenvalues of 
Ai for the Z^ group-cohomology wavefunction with the 
hexagonal branching structure, and I — 10,12,14. We 
find that the optimal fit for these eigenvalues is given by 
g = 0.925 ±0.01. 

One can also extract the central charges of these two 
CFTs by fitting the lowest A(/)i to —sq • I — ■ 1) 

|62) . In the Z 2 case we find c = 0.9983(1), and in the Z 3 
case we find c = 0.997(1). 

Lastly we comment on several 2-point and 4-point cor¬ 
relation functions that we obtained on a cylindrical ge¬ 
ometry. For the Z 3 case we find that charge neutrality 
has to be obeyed on the A sublattice and B + C sub¬ 
lattices separately, such that for example, the two point 
function between an 0+ operator on A and an 0_ op¬ 
erator on B is always zero. This is expected due to the 
extra Zjv symmetry implied by the hexagonal branching 
structure. Furthermore we find that (0+(r)0_(r')) is 
always positive when both are placed on the B sublat¬ 
tice (or when both are placed on the C sublattice) and 
always negative when one is on the B sublattice and the 
other on the C sublattice. Also considered are four point 
function ( 0 +(ra) 0 +(rf,) 0 _(r(,) 0 _(r(,)), on a cylindrical 
geometry, where by Vajh/c we denote both the position 
and the sublattice. Winding ri, around r^, we find that 
the complex phase of the correlator rotated by 27r as ex¬ 
pected. This behavior agrees very well with CFT predic¬ 
tions given the association between lattice operators and 
CFT operators made in the main text. 

For the Z 2 case, due to the extra Z 2 x Z 2 x Z 2 sym¬ 
metry, we find that charge neutrality has to be obeyed 
on each sublattice separately. In agreement with this 
fact and the two assumptions given in the main text, 
the following correspondence between microscopic and 
CFT operators is proposed: on the B (C) sublattice, 
O —>■ Ve=+i ± 14=-! (O —> i(t4=+i — Ve=-i)) and on 
the A sublattice, O —>■ Vm=+i + 14n=-i- This correspon¬ 
dence is again consistent with the previously mentioned 
numerical results. 


2. Non-Abelian 


We have performed the same transfer matrix numerical 
calculation for the case of The spectrum given in 
Table IV gives strong evidence for a free boson at g = 1. 




13 


CFT Operator 

A(3 = 1) 

A{1 = 16) 

A{1 = 18) 

A{1 = 20) 

I 

0 

0 

0 

0 

Vm=±l 

14=±i 

0.5(x4) 

0.50235616(x4) 

0.50186276(x4) 

0.50150942(x4) 

I {n = ±1) 

Ve = ihl,m = ±l 

l(x6) 

0.99901361(x6) 

0.9989935 (x6) 

0.99905296(x6) 

14i=±i <8) (n = ±1) 
14=±i (g) (n = ±1) 

1.5(x8) 

1.51518252(x8) 

1.51182344(x8) 

1.50947823(x8) 


TABLE IL Comparison between numerically obtained scaling dimensions and those of a compact boson with 
g = 1 , for the Z 2 group cohomology wavefunction. 


CFT Operator 

A(g = 0.925) 

A{1 = 10) 

A{1 = 12) 

A{1 = 14) 

I 

0 

0 

0 

0 

V^—±1 

0.4625(x2) 

0.46582886(x2) 

0.46457917(x2) 

0.46382596(x2) 

Ue = ±l 

0.540540(x2) 

0.54756391(x2) 

0.54568348(x2) 

0.5445524(x2) 

I ® {n = ±1) 

l(x2) 

1.0031008(x2) 

1.00042102(x2) 

0.99954883(x2) 

Ve=:±l,Tn=±l 

1.00304054(x4) 

1.00195587(x4) 

1.00156733(x4) 

1.0016539(x4) 

14,=±i ® {n = ±1) 

1.4625(x4) 

1.49370415(x4) 

1.48280295(x4) 

1.47672672(x4) 

14=±i (g) (n = ±1) 

1.54054054(x4) 

1.58217044(x4) 

1.56804518(x4) 

1.56029936(x4) 


TABLE III: Comparison between numerically obtained scaling dimensions and those of a compact boson with 
g = 0.925, for the Z3 group cohomology wavefunction. 


CFT Operator 

A(5 = 1) 

A{1 = 6 ) 

I 

0 

0 

14i=±i 

Ue=±l 

0.5(x4) 

0.518(0), 0.511(9), 0.513(0), 0.514(7) 

/ (g) (n = ±1) 
Ve=±l.m=il 

l(x6) 

0.949(7) (x2), 0.948(3) (x2), 0.949(0) (x2) 


TABLE IV: Comparison between numerically obtained scaling dimensions and those of a compact boson with 
g = 1, for the D 3 group cohomology wavefunction with a type 1 cocycle with p = 1 . 


Appendix B: Operator content and association for 
the c = — 7 loop model 

Here we analyze the operator content of the n = — 1 
dense loop model corresponding to the Levin-Gu wave- 
function. Our aim is both to identify which CFT opera¬ 
tors appear in this specific microscopic model as well as 
understand the association between this set of operators 
and the microscopic charge operators. To this end, we 
repeat the transfer matrix analysis used in the previous 
appendix, however this time on the full triangular lattice. 
The transfer matrix (T) generates the partition function 
(Z) for a cylinder of circumference I with one of the two 
primitive vectors of the lattice oriented along the circum¬ 
ference, such that straight lines of sites encircle it. The 
transfer matrix “propagates” along two such horizontal 
lines such that it does not mix the three sublattices. To 
allow a consistent labelling between the sublattices, I was 
taken to be a multiple of 3. 

The amplitude is simply given by (—1) to the 

number of domain walls in the conhguration. Since 


domain walls are not intersecting for Z 2 spins on a trian¬ 
gular lattice, they can be interpreted as loops and their 
number is well defined. The resulting 2D statistical me¬ 
chanics problem is a special case (n = —1, x = 1) of a 
general loop model whose partition function reads 

Z = (Bl) 

c 

where C is a loop configuration, ^[C] is the number of 
loops, l[C] is the total length of loops and x is the inverse 
loop tension. Under renormalization, it flows to a stable 
critical phase (the so-called dense-loop phase) with x = 

xo = (2- ^/2^)-l/2[5^1llHllSMl]■ 

This critical phase has a central chage of c = — 7 and 
is a non-unitary logarithmic CFT. Yet, based on the 
Coulomb gas approach, one expects primary operators 
to be some subset of electric charges (14) and 2Fleg wa¬ 
termelon operators {Wi) whose scaling dimension in this 
case are Ae = e(e — 2 eo)l 2 g and A; = — en/2ff, 

with Co = 2/3 and g = l/3 [i5l I3DH5^ . In Table [\^ we 
find a good match with the numerically obtained seal- 
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CFT Operator 

A(c = -7) 

A(/ = 18) 

K{1 = 21) 

CM 

II 

V^e=eo 

-2/3 

-2/3 

-2/3 

- 2/3 

II II 

^ A" 

-l/2(x2) 

-0.50090528(x2) 

-0.50067042(x2) 

-0.50051597(x2) 

III 

0(x2) 

0.00195229(x2) 

0.00141917(x2) 

0.00107904(x2) 

(dz,d4-Fe=eo 

l/3(x2) 

0.35635165(x2) 

0.35013343(x2) 

0.34614146(x2) 


TABLE V: Comparison between numerically obtained scaling dimensions for the Levin-Gu wavefunction and 
those of the Coulomb gas predictions for the c = —7 CFT. 


ing dimensions and the following subset of operators: 
Ve=eo,Wi=i,Ve=i,Ve=o = I,Wi =2 and the level one de¬ 
scendants of Ve^eo’ denoted {dz,dz)- 14 = 60 ■ These calcu¬ 
lations were performed for x = xq. 

Due to the presence of negative scaling dimensions, the 
previous method of extracting the central charge would 
now yield an effective central charge of c* = c — 12Amin. 
Under the above identification of operators, Amin = 
—2/3 and so we expect c* = —74-8 = 1. Numerically 
we obtain c* = 0.98(9). It should be noted that, due to 
the presence of negative scaling dimensions, one has to 
be careful about the boundary conditions at the ends of 
the cylinder. Coulomb gas results are only expected if 
one takes the identity CFT operator I as the boundary 
condition. 

Calculating {az{Q)<7z{r)) on long cylinders, we find an 
oscillation of the form cos(^r) (r = 0,1,2,...) times an 
exponential decay (for r > 1) with a decay length consis¬ 
tent with the scaling dimension of the operators Wi=i and 
14 = 1 - Assuming Wi=\ and 14=1 are related by the sub¬ 
lattice symmetry operations, the oscillation can be repro¬ 
duced by assuming the following correspondence between 
microscopic and CFT operators: 0{r) —)■ Oa = 14=i on 
the A sublattice, 0(r) —>■ Ob = (—1/2)14=1-I-•\/3/2W;=i 
on the B sublattice and 0(r) —> Oc = (— 1 / 2 ) 14=1 — 
on the C sublattice. 

There is obviously some arbitrariness in this choice of 
correspondence. This arbitrariness is related to the sub¬ 
lattice symmetries, as we now discuss, and has there¬ 
fore no physical impact, as it should. Let us define 
the 27r/3 rotation TZ. At the microscopic level, one has 
n{0{rA)) = 0(4^) where t'b = TZ{rA) and where the 
index indicates the sublattice to whom r belongs. At 
the CFT level, the counterpart of this relation would 
be TZ{OA{r)) = OB{r'). One can check from the corre¬ 
spondence above that this relation holds if one assumes 
that the two operators 14 =i and Wi=x form a real two- 
dimensional representation of the sublattice symmetry 
group such that 7^(14=i) = (—1/2)14=1 + (•\/3/2)Wi=i 
and 7^(W^=l) = (-v^/2)14=i - {l/2)Wi=i. Further as¬ 
suming that the sublattice symmetries commute with the 
conformal symmetry generators, one obtains that CFT 
correlators are invariant under the sublattice symmetries. 

Unlike in the ^Ar >2 group-cohomology case, the asso¬ 
ciation between lattice and CFT operators in the Levin- 


Gu case is such that lattice operators appear as a lin¬ 
ear combination of CFT primaries. The expression for 
a particular amplitude written in the charge basis would 
thus appear as a superposition of many CFT correlators. 
This complication is however avoidable following a simple 
change of basis. Pick one of the three main directions of 
the triangular lattice, say horizontal. For each horizontal 
BC bond, replace the local basis states \B) and \C) by 
|±) = ■^\B) ±\C). Leave the basis on the A sites as 
is. Now, the previous correspondence implies that 0(r) 
for |-|-) is proportional to I 4 0(r) for |—) is propor¬ 
tional to Wi while 0[r) for |A) is, as before, associated 
with 14- Using this basis, a charge configuration now 
translates into a single CFT correlator. 


Appendix C: The adiabatic connection with the 
dilute charge limit 

Here, evidence is given to show that the dilute charge 
limit, where non-zero charges are far apart on the lattice 
scale, is adiabatically connected to the dense charge limit. 
While the following approach is general, the numerical 
calculations were performed for the group cohomology [5] 
and the Levin-Gu [3] SPT wavefunctions. Our approach 
goes as follows. First, we consider the parent Hamiltoni¬ 
ans (Ho) of the original SPT wavefunctions |TH3], which 
contains no particle fugacity factor yet and therefore cor¬ 
respond to the dense charge limit. On closed surfaces, 
these Hamiltonians are local, gapped, and consist of a 
sum of commuting projectors (Hg = In Par¬ 

ticular, they are frustration-free Hamiltonians [66j : each 
term in the sum is positive definite and annihilates the 
ground state. Consequently, Hq can also be thought of as 
the parent Hamiltonians [66] associated with the Tensor 
Product State (TPS) version of these wavefunctions 
as shown in Eq. 

These TPS wavefunctions are then modified to 
by introducing a fugacity of j3 for any non-zero charge, 
such that /3 = 0 forbids charges altogether and /3 = 1 
yields the original wavefunction (|'0/3=i) = |'4o))- For¬ 
mally, this is done by using a non-unitary transforma¬ 
tion A^ = n./3'“ ’■I, where \ar\ is the operator measur¬ 
ing the (minimal) absolute value of the charge at site r 
(e.g. ar = N — 1 leads to \ar\ = | — 1| = 1) such that 
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\'4’p) = A^l'i/^o). For any finite /? this transformation has 
an inverse, and thus the following family of Hermitian 
local Hamiltonians can be defined: 

= (Cl) 



where is the product of on the sites r on which 
hi acts non-trivially. It can be easily verified that H,g 
also consists of local positive-definite operators which an¬ 
nihilate \'tpp) and thus is a parent Hamiltonian of IV'/?)- 
Except in the vicinity of /3 = 1, the Hamiltonians H^g are 
not guaranteed to be gapped. Nonetheless, provided that 
their ground states I'l/'/?) remain short-range correlated, it 
appears highly likely that they would be gapped. 

The correlation length of these modified TPSs can be 
studied analytically at the two extremities of /? = 0 and 
/3 = 1 and by diagonalizing the TPS transfer operator 
(see Fig. (§) at different circumferences (1) for gen¬ 
eral /3. Referring the reader to Ref. [67] as regards the 
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One simply obtains two copies of the previously used 
partition function (Z, Z), each with its own degrees of 
freedom and their corresponding “Boltzmann 

weights”, j). On top of these phase factors, 

the resulting doubled sum over “advanced” {(f)r) and “re¬ 
tarded” {(j)r) degrees of freedom now includes a factor of 
for each site, where = l—2(j)r. These terms 
represent the contraction of physical degrees of freedom 
between the retarded and advanced parts. Notably, at 
13 = 1 they force the retarded and advanced degrees of 
freedom to be locked together, while at /3 = 0 they leave 
them decoupled. The correlation length of the wavefunc- 
tion, with respect to any two local operators, is then 
bounded by the correlation length of this resulting “bi¬ 
layer” statistical mechanics problem. 


For /3 = 0 and /3 = 1, this correlation length (y) can be 
analytically evaluated. Note that for /3 = 1 the locking 
of retarded and advanced degrees of freedom removes all 
phases (since j = 1), and one obtains a sim¬ 

ple sum over only one set, say {(l>r}i without any phases. 
This is a random spin partition function with y = 0 as 
one expects from the original wavefunctions [2]. In the 
limit f3 ^ 0 the retarded and advanced parts become 
gradually decoupled and one obtains a doubled version 
of the CFTs studied in the main text. Consequently, % 
tends to infinity. For small but finite (3, a necessary (but 
not sufficient) condition to have a finite x is that charge 
fugacity be a relevant perturbation for the doubled-CFT. 
Scaling dimensions in the doubled-CFT are just twice 
those of the previously studied ones, and so we find that 
the scaling dimension of Uj 0 d'z is = —1 and 

^/ 3 ,Z 2 = 1 for the Levin-Gu and Z 2 -group-cohomology 
wavefunctions, respectively. In both cases it is smaller 
than 2 and consequently relevant as required. Based on 
this analysis, the simplest expectation for the phase di¬ 
agram is a single phase transition happening exactly at 
P — Q with a finite correlation length for any /3 > 0. 
Notably, the above argumentation generalizes to all Z^- 
group-cohomology wavefunctions. 


°0 0.01 0.02 0.03 0.04 

FIG. 5: Inverse correlation length (x~^) for the Z 2 group- 
cohomology (upper panel) and Levin-Gu wavefunctions (lower 
panel) on a cylinder of circumference I and for different val¬ 
ues of p. For Z 2 group-cohomology, the correlation length 
appears to remain finite, at least down to /3 = 0.2 where fi¬ 
nite size correlations begin to interfere. For the Levin-Gu 
wavefunction, x seems to decrease very rapidly with /3, which 
makes it possible to do finite-size calculations for smaller val¬ 
ues of P than in the group-cohomology case. 

technicalities of this procedure, the resulting numerical 
computation for the G = Z 2 cases is easily described. 


To corroborate the above expectations, we numerically 
obtained x~^ for cylinder circumferences of 1 = 6,9,12 by 
taking the difference of the logarithms of the two domi¬ 
nant eigenvalue of the TPS transfer operator. The results 
for the original Levin-Gu wavefunction (i.e. with x = 1) 
are shown in the lower panel of Fig. (§ where the pres¬ 
ence of a gap is evident even down to /3 = 0.05. For the 
Z 2 group cohomology wavefunctions shown in the upper 
panel, % appears to be generally larger thus limiting the 
regime were we can exclude a divergent x to /3 > 0.2. 
The X axis in Fig. ([^ was chosen to be one over the cir¬ 
cumference squared [3/P) since it gives the best linear 
fit. If we had chosen to plot the gap as a function of I/l 
instead, it would still appear that the gap remains finite 
above these values of p. 
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Appendix D: Corrections to the infinitely dilute 
charge limit 

The following expression was obtained for the 2^Ar>2 
group cohomology wavefunction: 

= Z-i f (Dl) 

i 

The action is given by 

5M = J + A(Fe=w + K=-^) 

+ y {Vm=N + Vrn.=-N) 

where only the most relevant vertex terms allowed by 
symmetry were included. In the main text, for simplic¬ 
ity, the infinitely dilute charge limit (p —>■ 0) was con¬ 
sidered, in which case A, A' —>■ 0 since these vertex terms 
are irrelevant. In that limit, since is non-zero 

only if i'ii® wavefunction possesses 

an enhanced {7(1) x t/(l) symmetry for the magnetic and 
electric charges instead of the microscopic Zn x sym¬ 
metry. 

As we show below, at small but finite p, A, A' be¬ 
come finite, and this acts to restore the original sym¬ 
metry. Furthermore, we argue that the dominant con¬ 
figurations are the ones for which the excess charges 
(Si 6* = Si i'O bunch in tightly bound 

A—particle clouds whose core size is of the lattice scale. 
Under this picture, the wavefunction retains the factor- 
izable form appearing in the main text and one can es¬ 
tablish hidden order at all scales. 

Let I denote the average distance between charges, such 
that p = UqIP. Between the lattice scale (oq) and I, one 
may apply the usual RG procedure on the action, and 
find that A and A' are suppressed as , with 

A being N'^/2g and gN‘^/2, respectively. For neutral 
configurations {E,M = 0), when charges are far apart 
on the lattice scale, the finiteness of A typically plays no 
role. The strong attraction between +N and —N charges 
would make them form neutral pairs on a microscopic 
scale [65 such that they would not interfere with the 
long range behavior. 

For a non-neutral configuration with E = N, X be¬ 
comes essential for obeying charge neutrality in the 
CFT correlator, thus making the corresponding ampli¬ 
tude non-zero. This non-zero probability weight in the 
E = kN sectors with k G Z\0 restores the original 
symmetry associated with electric particles. In the fol¬ 
lowing, we discuss an electric charge imbalance for sim¬ 
plicity, but the exact same argument applies to a mag¬ 
netic charge imbalance. 

For notational purposes, it is advantageous to enumer¬ 
ate the set {ei,mi,ri} by the positions of the N+e 
particles of charge e = -|-I, the positions zj of the N_e 
particles of charge e = — 1, the positions w'^ of the N+m 
particles of charge m = -l-I and the positions w'f of the 


N-m particles of charge m = — 1. With this notation, 
one has -z/; Wj!"; rC;”), E = 

and M = Nj^^n — N_m- We now look at the wavefunction 
amplitudes in the E = N sector. At first order in A, they 
are given by 



{Ve^l{z+)Ve=-l{Zj )Vm=l{w+)Vm=-l{Wi )Ue=-Ar(s))o 

(D3) 

where ■ ■)q corresponds to taking the average for A, A' = 
0. While the integrand can still be expressed in a 
Laughlin-like product form, the integration over the po¬ 
sition of the charge —N operator (henceforth referred to 
as the dust particle) makes the total amplitude unfactor- 
izable. 

Notwithstanding, we conjecture that typically the 
above amplitude is in fact factorizable and retains a sim¬ 
ple Laughlin-like form. As argued below, typical configu¬ 
rations are ones in which the dust particle is surrounded 
by a tight cloud made out of the excess E = N charges. 
The core of this cloud is of the size of the microscopic 
lattice scale (oq) and the integration over the position of 
the dust is effectively limited to this microscopic region 
and may be removed. Furthermore, the core of the cloud 
renormalizes the charge of the dust particle down to some 
effective value Neff. As a result, the Laughlin-like form 
for the wavefunction is restored in a “renormalized” form 
with new degrees of freedom consisting of the renormal¬ 
ized dust particle and all the original physical particles, 
excluding those in the cloud’s core. Consequently, the 
arguments for hidden order used in the main text can be 
again applied using the plasma derived from this renor¬ 
malized wavefunction. 

To support the above conjecture, we wish to show that 
it is self-consistent in terms of what the dominant config¬ 
urations are. Dominant configurations correspond to a 
large magnitude \A{zl' ; z~ ; w'^ ; wf ) |. Let us assume this 
magnitude can be written as 



\{Ve^l{zt)Ve^-l{z^ )U„=l(w)!')Um=-l(W/ )Ue=_jv(s))| 

(D4) 

which is allowed provided that the integral over (Ps 
is dominated by a set of configurations for which the 
complex phase of A{zN, z~ ] ; wf ) varies very little, 
as assumed previously. We will then show that un¬ 
der this assumption, the dominant configurations of 
\A{zf ; z~ ; w '^; rc;~) | are those in which it is indeed of this 
form. 

The main advantage of Eq. |D4| is that the integrand in 
the right hand side now appears as the Boltzmann weight 
of two decoupled plasmas associated with the magnetic 
and electric particles. Relevant to our discussion is the 
electric plasma which, in plasma terminology [6 1 j . is a 
two-component plasma with T = 1/g and one dust parti- 
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cle of charge —N. For such values of F <2, the plasma is 
screening with a screening length of the order of 1. Conse¬ 
quently, dominant configurations are those in which the 
excess charge form a screening cloud around the dust 
particle. 

Let us estimate the form of this screening cloud for 
large N. To start with, we consider the subsystem con¬ 
sisting of the (negative) dust particle and a single positive 
charge. The partition function of this subsystem is 


■^dust+l 


drs 


(Pzi 

al \s-Zi\^/3' 


(D5) 


For N/g > 2, the ultraviolet divergence in the above par¬ 
tition function implies that the partition function is dom¬ 
inated by configurations where zi is bound to s within 
a small microscopic region (ap). Effectively, these par¬ 
ticles are then paired into a composite object with a 
charge of —N + 1. We may repeat the argument with 
another charge 1 particle which, provided the condition 
{N — l)/g > 2 is satisfied, will bind itself to the com¬ 
posite object. Consequently, already at the microscopic 
scale, the dust charge will be strongly reduced down to 
—N^g = — \2g\ since it would have N — N^g charge -1-1 
particles bound to it. 


We may now easily justify assumption D4 self- 
consistently. A configuration {zf ; z~]w '^; wf) is deemed 
typical if — A^eff of the A^+e particles present at posi¬ 
tions z^ form a tightly bound cluster. It is then advanta¬ 
geous to separate the set of charge -1-1 particles that form 
this cluster {z^^ with p = 1,... ,N — Neg) from the other 
charge -1-1 particles with q = 1,..., N+e — A^ -I- A^eff)- 
With this notation, one can write \z^^ — z\ zz ao W p 
where z is the center of mass of the positions z^^, i.e. 
z = ~ ^es)- For a typical configuration, the 

integral over d?s is only appreciable within the bunch¬ 
ing region, i.e. within a distance of the order of oq 
from z. Since magnetic charges appear on the scale of 
Z ^ oo, typically they would not enter this small re¬ 
gion. Consequently the variation of the complex phase 
of the integrand scales as oq /I within the effective region 


of integration. Thus, in the limit of small Og/Z, Eq. D4 


is justified. Finally, the effective Laughlin-like form for 
A{z^^,z^^;z~;w'j^;w'f) in the E = N sector is simply 
given by 


M4p > 4 ; ) = ipiz^ - 5) H ^ 

iqjkl 

(Ve=i(z+)Ve=_i(z- )Vm=i{w+ ) 14=-^,//(z)) 

(D6) 

where ip{z^ — z) is the short-scale wavefunction for the 
particles in the cluster and has no impact on the long- 
range behavior. Now that the wavefunction amplitudes 
are factorized, it is possible to show hidden order in the 
same way as done in the main text. 


Appendix E: Relation to previous work 

As explained in the main introduction, some field the¬ 
ory arguments were given in order to identify the edge 
theories of SPT states for which a continuous space-time 
description applies. Two different techniques were used. 
The first one relies on embedding the symmetry group G 
in (a product of several of) Z7(l) group(s) and to study 
the associated Chern-Simons theories with unimodular 
K matrices|llL IM|. This technique is limited to 2-l-lD 
dimensions. Abelian groups and to cocycles that have a 
simple Chern-Simons interpretation (so-called type 1 and 
type 2 cocycles) [53] ■ For all the cases treated in these ref¬ 
erences, the obtained edge theory was a non-chiral free 
boson (c = 1) with symmetry breaking terms. Since the 
relevance or irrelevance of these terms depends on the 
Luttinger parameter, which is non-universal and depends 
on the microscopies, this argument basically states the 
edge theory is either a free boson or a symmetry-broken 
phase. Interestingly, our microscopic calculations show 
that, for the cases we have looked at, the Luttinger pa¬ 
rameter is in a regime where the theory is critical, which 
was not guaranteed. 

The second one relies on embedding the symmetry 
group G in SO{D+2) (for a SPT in D spatial dimensions) 
and to look at a non-linear sigma model with 9 — 2'Kk 
with fc € Z [TS] |3D] . This technique is limited to groups 
that can be embedded in SO{D 2). There, for 2-1-1- 
dimensional SPTs, the predicted edge theory is SU{2)i 
plus some symmetry breaking terms. 

We want to emphasize that, unlike these field theory 
techniques, (1) our technique applies deep in the SPT 
phase and (2) it applies to groups too big to be embedded 
in SO{D + 2) (like Z| for D = 2), to non-Abelian groups 
and to type 3 cocycles. We have at least one example 
{G = Zg with a type 3 cocycle) for which c = 2 and for 
which it is not obvious how to deduce it from these field 
theory arguments. 


Appendix F: Branching structnre 

Quite surprisingly, although edges of SPT states can be 
gapped (if they are symmetry-broken), all the entangle¬ 
ment spectra obtained here were gapless. Furthermore, 
as discussed earlier, some are even integrable. In general, 
one may wonder whether this is a universal property of 
such wave functions or, equivalently, whether the strange 
correlator of an SPT state with respect to a trivial state 
is always critical. While we still lack a comprehensive 
answer to this question, it seems that in all cases an¬ 
alyzed by us and others [30], critical behavior was ob¬ 
tained. Notwithstanding, some changes to the wavefunc¬ 
tion do result in changes to its corresponding CFT or at 
least to its convergence properties. 

One example of such a change is a chance in the 
branching structure. Instead of the ABC branching 
structure used throughout this work, let us consider us- 
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ing a simpler branching structure for which the repeat¬ 
ing unit is a parallelogram spanning four sites. It can 
be shown that for the case oi G = Z 2 and the choice of 
cocycle we have made, this has no effect. However for 
G = Z 3 , it does alter the partition function. Repeating 
the numerics for such a branching structure still results 
in a critical theory (i.e. transfer matrix eigenvalues scale 
roughly as 1 /L) but the scaling dimensions seem far from 
converging to their asymptotic CFT values. For instance, 
descendants were not clearly distinguishable in the spec¬ 
trum for the system sizes we could reach. 

If one were to alter the branching structure very 
sparsely, this behavior could be understood from the 
CFT perspective developed here. The ABC branching 
structure resulted in an extra Zj^ symmetry acting only 
on the A sublattice. This extra symmetry is removed 
by the altered branching structure. This allows terms 
such as cos{ 4 > — 0 ) in the action having a scaling dimen¬ 
sion of l/2g + g 12. These terms are relevant and should 
therefore drive the model away from the c = 1 theory, to 
either a different critical theory or a gapped, symmetry- 
broken theory [70]. Notably, for Z 3 , we have verified that 
the altered branching structure still results in a critical 
theory. 

Thus the ABC branching structure which we have used 
throughout this text, and probably also the cocycle itself, 
has an important consequence on both the analytical and 
numerical tractability of the resulting CFTs. While ex¬ 
pected from the technical perspective, this is quite un¬ 
usual from the mathematical perspective where this free¬ 
dom of choice should have no physical consequences and 
in which there are no canonical choices. 


Appendix G: Flux responses 

In an Integer Quantum Hall phase with a^y = e^//i, 
inserting half a flux quantum draws in half of an electron. 
Considering SPT phases based on type ii cocycles, simi¬ 
lar quantized responses to symmetry fluxes have been re¬ 
cently identified m- Here we show how this is reflected 
in our formalism. 

To study such effects, static symmetry fluxes should be 
introduced into the formalism. To this end, we consider 
a disk geometry, we add a (static) gauge field A; G G on 
each link I of lattice and we replace in the Hamiltonian 
all terms such as (pi — (pj by (pi — pj + Ai ^., where Uj is 
the link connecting site i to site j. Without fluxes. A; 
can be written as the discrete derivative of a scalar fi 
(i.e. A;.^. = fj — fi) and can thus be removed using a 
local gauge transformation taking pi ^ pi — fi on each 
site. A flux g G G at the centre of the disk means that 
Ai = g where O is any set of links surrounding the 
origin. Note also the trivial fact that every path that 
does not encircle the origin still obeys ~ 

Next we wish to carry flux insertion from the Hamil¬ 
tonian level to the ground state wavefunction level, or in 
other words relate the ground state wavefunctions before 




FIG. 6: Introducing static symmetry fluxes into the wave 
function. Acting with a unitary transformation which rotates 
by g the right hand side of the system can be thought of as 
two opposite fluxes at the origin (see (a) and (&)). Removing 
the lower flux yields the g-twisted wave function (see (c)). 
The location of the branch cut is insignificant and may be 
shifted using a unitary transformation (see (d)). 


and after the flux insertion. Deferring proofs for a mo¬ 
ment, the prescription for doing so is the following: First 
consider both g and g~^ fluxes, or equivalently act with a 
unitary transformation which rotates by g all sites which 
are, say, on the right part o f the disk (see Fig. (l 6 l)(a,b)). 
As shown in Section IV A (see also Ref. EH), the phys¬ 
ical symmetry simply rotates the flelds of the CFT {tp 
and its dual 9) in a certain cocycle-dependent way. Thus 
the CFT obtained following this unitary transformation 
would simply have a 5 -twist in the boundary conditions 
of the CFT flelds when crossing the upper vertical and 
5 “^-twist when crossing the lower line. Since the effects 
of the 5 -twist and of the 5 “^-twist are distinct, one can 
now simply omit the 5 “^-twist (see Fig. (|^(c)). The 
resulting wave function is what we call the 5 —twisted 
groundstate wavefunction. 

Next we explain how the 5 -twisted wave function is re¬ 
lated to the groundstate of the Hamiltonian in the pres¬ 
ence of a 5 flux. Let us discuss the case which is of 
relevance here where the Hamiltonian is a sum of local 
projectors which annihilate the ground state. Introduc¬ 
ing a 5 flux in the Hamiltonian as in Fig. Q(c) amounts 
to altering the projectors which act on both sides of this 
line via a Peierls-like substitution. Other projectors re¬ 
main unaltered and thus a few correlation lengths (y) 
away from the 5 —flux insertion line, projectors acting on 
the 5 -twisted ground state still annihilate it. The reason 
is that, in a gapped system, all effects are local and thus 
the density matrix of the 5 —twisted groundstate in this 
region is exponentially close to the original ground state. 

Next we consider the subset of projectors which are 
close to the flux insertion line (on the scale of y) but 
still far from the origin. Here it is beneficial to use a 
unitary transformation and shift the flux line away from 
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these projectors (see Fig (|^(d)). This removes the flux 
from both the projectors themselves and also, due to the 
finite correlation length, from the wave function. Af¬ 
ter this transformation, it is clear that the current sub¬ 
set of projectors again annihilates the g-twisted ground- 
state. However, since this transformation is unitary, it 
also means that, even prior to acting with it, the subset 
of projectors annihilated the g—twisted ground-state. Of 
course, it must then be that, also prior to this unitary 
transformation, it was annihilated by these projectors. 
We thus find that all projectors annihilate the g—twisted 
groundstate except perhaps those which are a few corre¬ 
lation lengths away from the flux insertion point. This 
implies that the g-twisted groundstate wave function cap¬ 
tures the long range effect of the flux insertion (see also 
Ref. EH). 


Consider G = with the simplest type ii co¬ 

cycle (pi = pii = 0 ,p' = 1 ) and pick g which acts only 
on the first group (g = gi x I). As shown in Section 
IV A within our CFT formulation the physical symme¬ 


try transformation g rotates only the 9 degree of freedom 
of the compact boson by The 5 -twisted wavefunc- 
tion thus has a CFT with a g—twist only for the 9 field. 
A few correlation lengths away from the flux insertion 
points, the effect of the flux thus appears as an insertion 
of 1 /A^ of the basic magnetic charge. 

In the dilute charge limit, the effect of inserting a frac¬ 
tional magnetic charge can be deduced from the plasma 
analogy used in the main text. From the two two- 
component plasmas associated with the norm of the wave 
function, the electric one would remain unchanged. How¬ 
ever the magnetic one would contain an external point 
charge at the origin with a fraction of the fundamental 
charge. This external charge would get screened by the 
physical charges. When measuring the charge within a 
large region around the origin (or more accurately the 
difference in this quantity before and after the flux in¬ 
sertion), the external point charge would not be counted 


and only its screening cloud would be picked up, yielding 
a fractional result. 


Considering a generic SPT phase with G = Z^ and a 
type i cocycle, a complication arises. As shown in Section 


IV A the physical symmetry now rotates both ip and 9. 


Within the plasma picture, this amounts to adding both 
half a magnetic and half an electric charge at the origin. 
Both electric and magnetic plasmas would then screen 
the external charge relevant to them and two screening 
clouds would be formed: One carrying half an electric 
charge and one carrying half a magnetic charge. Since 
the physical symmetry measures their sum, no fractional 
response associated with G = Zm would be detected. In¬ 
deed, in Ref. (EH only type ii cocycles were associated 
with fractional responses. Its worth mentioning that the 
specific type i wave functions considered in this work, 
which enjoy an enhanced Zj^ x Zn symmetry, would still 
see fractional charges associated with this enhanced sym¬ 
metry. In particular, measuring the total charge of the 
Zn symmetry acting on the A sublattice would pick up 
only the magnetic screening cloud and thus show a sim¬ 
ilar fractionalization effect as before. 


Based on the above result, we conjecture that a resid¬ 
ual quantization effect remains even if the enhanced 
Zjv X Zjv symmetry is broken down to only the diag¬ 
onal Zn- Let niy be the operator which counts symme¬ 
try charges associated with only one of the ^at’s on site 
i. Let Ni^i be the sum of niy within a radius I 3> x 
around the flux insertion point. As previously argued, 
the difference of Nij before and after the flux insertion 
{ANi,i) would be fractional for the enhanced symmetry. 
Next consider a perturbation V breaking Z]\j x Zjq down 
to its diagonal part and present the wave function as a 
superposition of A^i^oo eigenvalues. Then consider the op¬ 
erator Ni^i mod 1 . Importantly, perturbation theory 
in V only generates corrections to Ni^i that are integer, 
and therefore A{Ni i mod 1 ) would still be fractional. 
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